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ABSTRACT: We find a Umit cycle in a quasi-equilibrium model of evaporative cooling of a 
two-component trapped fermion gas. The existence of such a limit cycle represents an obstruction 
to reaching the quantum ground state evaporatively. We show that evaporatively ~ 0{1). We 
speculate that one may be able to cool an atomic fermi gas further by photoasso elating dimers near 
the bottom of the fermi sea. 



I. Introduction: The spectacular successes of laser cooling techniques in creating Bose-Einstein 
condensation (BEC) in trapped dilute alkali vapors'^'^'^l has stimulated efforts to form dilute nearly- 
degenerate atomic fermion gases. Such systems undoubtedly have unique phenomenology, and since 
one can control the composition, densities and even scattering lengths in principle, they furnish a win- 
dow to familiar phenomena (such as superconductivity, etc.) in unusual parameter regimes. t^'^'^'^'^l 
A critical step in achieving BEC in dilute alkali vapors is evaporative cooling. In this note we de- 
scribe limitations to the use of evaporative cooling for a harmonically trapped two-component fermi 
system. 

Overall antisymmetry of the final state wavefunction forces the s-wave scattering amplitude for 
two-body collisions in a single component polarized fermi gas to vanish identically. However, in a 
system composed of two or more fermi species there can still be appreciable s-wave scattering ampli- 
tudes at low energies. Recent experimental observation of quantum statistical effects reducing the 
scattering frequency at low temperature has been reported in reference [9] . Sympathetic evaporative 
cooling in two-component Bose systems has been experimentally verified^°. Aspects of the dynamics 
of sympathetic cooling in a two component fermi system have been discussed theoretically^^'^^, and 
recently achieved experimentally-^^. 

We will show in a robust model that sympathetic cooling of trapped fermions is intrinsically 
limited to ^ ~ where T is the temperature and fi is the chemical potential. Largely independent 
of trap and atomic parameters, this limit indicates that sympathetic evaporative cooling alone cannot 
achieve occupation probabilities in the trap single particle ground state characteristic of typical 
degenerate fermi systems (e.g. atomic nuclei and typical metals). 

In summary we model evaporation as a succession of quasi-equilibrium states. We find that 
evaporation moves the chemical potential towards saturation at a fixed (non-zero) fraction of the 
evaporation energy. Lowering the evaporation energy in an attempt to cool further simply causes 
the fermi surface and temperature to both recede, thus not substantially increasing the occupation 
probability of the lowest single particle state. 

We demonstrate this limiting behavior in a model with two regimes, one for which the evapora- 
tive energy scales of both fermion types are identical and the opposite extreme where one selectively 
only allows evaporation of one species from the trap. Rather than end on a pessimistic note, we 
conclude by speculating on a possible method for surmounting the difficulty of cooling an atomic 
fermi gas. 

II. The Model: We aim to explain general features of the evaporative cooling of a two-component 
fermion gas without recourse to the details of dynamics and transport. In particular, trap lifetime 
and other timescales will play almost no role in our considerations. 

As described above, we focus on the case of a two-fermion system cooling sympathetically by 
evaporation through interspecial two-body collisions only. As described in the introduction, this is 
a reasonable assumption as long as the density is low and the temperature is low enough that the 
amplitude for intraspecial scattering is small. 

In the limit of interest for the first regime, we model the e-vaporation process as a succession 
of quasi-equilibrium states in which the distribution functions of both species are cut off at Egyap- 
Aspects of the thermodynamics of dilute harmonically trapped fermions without an energy cut- 
off are described in Ref.[5]. Furthermore, we ignore the effect of any other environmental fields 



(for example, trap fields), and assume them to be constant over the lifetime of the system. The 
average effect of the interactions between the different species is absorbed into a mean field term 
which we sweep into fi (see, for example, Ref. [18]). We further assume that particle number is 
not communicated between the species, and so their individual /x can differ. Interspecial scattering 
processes do communicate energy between species, and so we assume that both species are always a 
at common temperature 

Consider evaporating both species at the same energy cut-off E^vap- We complete the analysis 
for this case and then turn to the opposite extreme where only one species evaporates. If both 
fermion species a and h are (nearly) the same mass then crossing symmetry equates the total rate 
for scattering into final state \a,h> with that of \h,a >. This means that for the same evaporative 
cutoff, the rate at which particles of one species evaporates equals that of the other. 

Finally, for simplicity, we model the evaporative process as one which always reduces the particle 
number by one and removes energy E^^ap- In the model we develop we will ignore the contribution 
to the cooling that results from the interspecial mean field^^. In the cases of interest in current 
experiments this interspecial mean field energy is expected to be very small compared with the 
other relevant energy scales (for example the fermi energy). Although this model of evaporation is a 
gross simplification, it becomes a progressively better approximation as the temperature drops, and 
we are confident it captures the main features of the evaporative process. 

We approximate each component's scaled number and energy by cut-off equilibrium distribution 
functions 

~ Jo e/3(^-M) + 1 ~ Jq e^(^-^^) + 1 ^ '* 

where (3 and are respectively the inverse temperature and the chemical potential both made 
dimensionless by factors of E^vap- Further, d depends on the actual spectrum of the trap, and for 
a three dimensional isotropic harmonic trap is 2. We will keep the discussion rather general with 
respect to d, but use d = 2 in all the graphs and particular conclusions below. While the actual 
value of d does not substantially impact the nature of the conclusions that we draw, we do find that 
reducing d (by, for example, significantly changing the aspect ratio of the trap) makes sympathetic 
evaporative cooling generally less effective. Finally in Eq.(l) E is dimensionless {i.e. in units of 
Eevap) and we have suppressed some overall factors that depend on E^^ap £^nd trap frequency. 
We model evaporative cooling by simply following Eq.(l) through flow along 

By using symmetries of the scattering for fermion species of nearly equal mass for the case where 
both fermion species are being evaporated at the same energy E^yap, we find that the net effect on 
the individual distribution functions is encapsulated in Eq.(l) and Eq.(2) for each species separately 
(and so in what follows for this regime we suppress indices). 

The resulting differential equations for (3 and n along this evaporative trajectory read, 
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where the determinant det{M) is given via 



(5) 



which, from the quasi-equihbrium distribution functions Eq.(l) we find det{M) > for aU d and 
This positivity may be understood on general grounds via the connection between det{M) and 
the specific heat, cy at constant N, 

P^detiM) 

= dN] (6) 

and by the fact that > as a consequence of, for example, Eq.(l). The det{M) vanishes in 

2 4 

the low temperature limit as ~ See Figure 1 for an example of this behavior (for d = 2 and 

/x = 2/3). 

log,„det[M] -vs- 



Figure 1: det{M) for fi = 2/3, as a function of p. 

Starting far from degeneracy, the net effect of the evaporative process is to alter fj. while in- 
creasing (3. Note that using the equilibrium distributions Eq.(l) implies by Eq.(3) that (3 increases 
monotonically along the flow. Upon reflection we see that, //'s evolution does not share this property. 
Instead, we flnd that the system Eqs.(l), (3) and (4) has a limit cycle at ^|^ — ^|^ = Oin the 
(/?,//) plane. This limit cycle is at an intermediate value of n — for all temperatures. 

Figure 2 below is a graph of //* as a function of temperature for d = 2 

Limiting Critlcol Chemlcol Potentlol Versus Temperoture 
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Figure 2: The Evaporative Limit Cycle with the Low- Temperature Asymptote 

For d = 2, the low temperature limit of fj,* is 2/3. We reiterate, that for values of 
below the curve, the evaporative process increases whereas for (/?, //) values above the curve they 
reduce it. Thus, as one tries to lower the temperature lower by reducing Eevap (thus increasing fx), 



evaporation causes the fermi surface ii*Ef,vap to also drop away. Of course, as one drops E^vap the 
scaled inverse temperature /3 also drops trivially. 

A figure of merit measuring how close one is to the quantum many body ground state is the 
occupation probability in the trap single particle ground state. This is a function of the ratio of 
the chemical potential to the temperature, in our notation Pfi. This combination is independent 
of Eevap in our simple model. We now show that this product is limited by the total scaled atom 
number (which in the normalization of Eq.(l) is limited to be at most -j^)- 

We consider two cooling methodologies which we refer to alternatively as "passive" and "ac- 
tive". They refer respectively to holding the E^vap fixed or suddenly reducing it. Unfortunately, 
in actual experiments underway E^vap is varied continuously, so these two cooling methodologies 
are probably not directly related to current experiments. This choice of cooling methodologies is 
advantageous analytically since it does allow evaluation of the effects of evaporation without recourse 
to any dynamical timescales. 

In the passive method, the system is held at a fixed E^vap and allowed to cool by evaporation 
indefinitely. Dynamically, since the emerging fermi surface at //* is always a fraction of E^yap, the 
cooling rate is limited by the overall rate of escape, which (at low temperature) is Boltzmann sup- 
pressed by a factor of e~^'^^~^*\ and by the trap lifetime. However, we reiterate the "kinetic" model 
we employ analytically encapsulates limits to sympathetic evaporative cooling of the two component 
fermion system without including "dynamical" effects (such as coUisional or trap timescales) . 

In the active method on the other hand we abruptly lower Eevap to a value at or below our 
initial iJ,*Eevap- Call this new evaporation energy E'^^^p. The distribution function is invariant under 
such a change in Eg^ap- The fj, and P values will trivially jump by factors of the evaporative energy 
scale ratio. Of course, the overall scaled phase space constants that we suppressed in Eq.(l) do go as 
positive powers of E^vap- Dropping Eevap to E'^^^p < Eevap has the immediate effect of dropping N, 
for example. This corresponds precisely to the statement that the all the particles with E > E'^^^p 
leave immediately. 

Since nfi measures our progress towards the ground state and is invariant under a sudden drop 
in Eevap-, we see that only subsequent evaporation (and rethermalization) of the remaining fermions 
can result in increases in the product We now show that evaporation after the drop in Eevap 
does not lead the system substantially closer to the quantum ground state. 

To recapitulate, the relevant question in both cooling methodologies considered here is whether 
the subsequent evaporative cooling leads to a large product. 

The evaporation equations above can be integrated numerically for /?(/i) by eliminating N 
from Eq.(3) and Eq.(4). In doing so, one finds that for values of ^ /j,*, the change in the 
dimensionless temperature ratio /? is generally relatively small, on the order of /3 itself. To get 
to a nearly degenerate fermi system starting far away from the ground state, we need a cooling 
regime in which much larger temperature drops are achievable. In numerical simulations, one finds 
that the only large temperature changes happen evaporatively when the system is at a |U very near 
//*, basically within 10 percent of that value. Figure 3 is the integral (3{n) of Eq.(3) and Eq.(4) 
for two initial conditions, {iii, j3i) = (0,2.5) and (2.8,4.5). Trajectories that start at higher initial 
temperatures (lower Pi) remain substantially lower on Figure 3 throughout the entire evaporative 
trajectory, but do eventually wind along the limit cycle at n*. towards large p. Recall that in the 
active method we are actually starting generally at smaller Pi than shown. 
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Figure 3: Two Typical Cooling Trajectories, both starting at P 4, from above and below /i* . 

Curiously, note that for /j, » the product can actually initially decrease along the 
evaporative trajectory. In the conclusion we comment on a heuristic way of understanding such 
counter-intuitive behavior. At any rate, it is clear from Figure 3 that achieving a condensed fermi 
system requires evolving /x close to //*. We now investigate some liberal bounds on how well evapo- 
ration can achieve that goal. 

In Figure 3 we integrated Eq.(3) and Eq.(4) for by eliminating N, not taking into account 
that N must always be decreasing! We now study how far along these trajectories we may 
progress evaporatively until we substantially run out of particles. 

Analytically, at low temperature, we find that the fiow equation for /j, near //* is independent 
of P entirely, and reads. 
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which can be integrated for all /i to read 
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The number evaporated AN must of course be less than the total number of particles in the trap. 
Using Eq.(l), we see that there are indeed stringent limits on the R.H.S. of Eq.(8). Since we 
know that appreciable cooling in this scheme doesn't occur until one is close to /U*(=2/3 at low 
temperature), we know that the best one can do is to evaporate all the particles in excess of the 
ground state at //*. For the scenario in which we start at a above //* (for example, as may be 
created in the active method) this limits the RHS substantially. For /i < //*, the system evaporatively 
evolves towards hut can never reach because there are simply not enough particles to evaporate. 
That is one reason to use the active method in a phase of the cooling, since it can raise the initial 
ji above the ji*. However, raising by this means is also self- limiting for two reasons. First, if 
it is raised substantially above 1 too many particles are lost from the trap and there are too few 
remaining to evaporate back to degeneracy. Secondly, as described earlier, it reduces Pi by the same 
factor it increases /U^, indicating the need to get even closer to to wind along the limit cycle and 
recover large p. We explain this in more detail quantitatively below. 

We graphically describe the consequences of Eq.(8). All the discussion here is in the "best 
case" scenario, in that we imagine starting with a system at already relatively large P (so Eq.(8) 
applies) and ask how well evaporative cooling can further increase P and thus Pii. Figure 4 is a 
graph of Eq.(8). The dashed line represents a maximum possible RHS at that initial ji. As per the 



preceding discussion we have plotted the contribution from the total number of particles for 
and plot the excess only for /i > /i*. The light dotted trace is the LHS of Eq.(8). Thus, to estimate 
the maximum possible increase/decrease in |U, start at some initial use the height of the dashed 
line to estimate how much of a change in the height of the light dotted line you may achieve. The 
resulting position at that height on the light dotted curve then is an (over-)estimate of the largest 
jjL achievable! 

Graph of Eq.(8) 



Figure 4-' Graph of Eq.(8) 

Studying the graph indicates that the process allows one to get within perhaps 10% of //* at 
best. As a very crude estimate, we can see that this occasions at most a roughly 3-fold increase in (3. 
Cooling from a typical active-method initial state of (/?, jj) ~ (1, 1), we find that one cannot achieve 
products in excess of roughly 3. For such a gas, the typical occupation probability of the lowest 
energy 1-particle state in the trap is roughly 90%. Of course, measuring this occupation probability 
is not a good measure of quantum degeneracy because it is exponentially sensitive to We prefer 
thinking about nearly degenerate quantum fermi systems as those with large products such as 
are found for nucleons in nuclei and electrons in typical metallic systems. 



Ill: A Second Model: We now study the scenario in which the evaporative thresholds for the two 
species are very different. We forgo a detailed quantitative analysis specific for this case, and instead 
reduce to and reason from the simpler model we have already discussed. For simplicity we assume 
that the masses and the trap potential of the two species are identical. Thus, the quasi-equilibrium 
expectation values of Na,b and Ea,b are given by the obvious doubling- and indexing- of Eq.(l). 
However to include the fact that the evaporative threshold for species a to be so much larger than 
that of species h replace the upper limit '1' of the integrals in the Eq.(l) by 'oo'. We then scale /3, 
and lib by the evaporative threshold of the h species. Subsequently E^vap refers to the evaporative 
threshold for species b only. 

The condition, dA^a = now becomes a linear constraint in the space {/la, fJ-b, P)- All these 
variables are dimensionless, scaled by the appropriate factors of E^vap- 

That linear constraint reduces the evaporative evolution of this two-component fermi system 
again to a two-dimensional dynamical system {i.e. at fixed E^vap and Na, is really a function 
of (3). In particular, scattering events that leads to evaporation of a particle of species h at energy 
Eevap will, on average, remove a net Ea from the total energy of species a and E^ from species b, 

t Note that although is always greater than zero, the total derivative ^| < for // < //* 

because ^l/j,^ is negative. 



where Ea + Ei, = E^yap ^ind with ratio Ea/E^, depending on /3). In equations, dAT^ = 

imphes that 

dEa ^ detjMq) 

(compare to Eq.(6)) where det{Ma) is the determinant of the matrix of partial derivatives for the a 
system only. Note that this determinant det{Ma) is now computed with the integrals extending to 
oo, and so is singular at high temperatures, but still looks like the rest of the graph in Figure 1 for 
low temperatures. The analogous function for the h species, det{Mi,), is precisely the same as for 
Eq.(5) with Eq.(l) (that is, using integration limits [0 : 1]). 

Energy conservation implies (in scaled dimensionless quantities) that the evaporative trajectory 
is along 

div, (10) 
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The differential relations between Eh, Nh, lib and /3 are exactly the same as for the matrix system 
analyzed in the first model. We now use Eq.(9) to rewrite the dE^ term in Eq.(lO) in terms of dj3 
and use it to re- write the two dimensional system for the evolution of and (3 in terms of dA^^. 
Recall that Ha also changes, but is given parametrically in terms of (3 (and Na-, which is held fixed). 
We find, 
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The denominator is 



det{M) = det{Mb) + det{Ma) (13) 



and is again strictly positive. 

We can now use the intuition based on the first model discussed in detail above to constrain 
how well one can cool the a species in this scenario. The measure of how close to the ground state 
we are for that species is again the product f3^a- Consider two possible initial high temperature 
limits: Na » Nb and Na « Nb- In the first case, Ha is high but there are fewer b species to 
evaporate off, and by the analysis of the first model, we expect only a modest increase in /3. In 
the other extreme, /Xa starts off small, and so, although in principle you could carry off much heat, 
it's temperature will track the temperature of the b species, which will again be limited by all the 
considerations discussed in the first model, and so the product will start small and stay small. 
Thus, the most promising initial (high temperature) state is Na ~ Nb (that is, ~ /^b). 

Looking now at Eq.(12), we see that the limit cycle still exists and occurs at a value //^ lower 
than in the first model we discussed. Also note that the /? evolution equation in this scenario 
differs only by the prefactor, det{Ai), from the first model we considered. Due to positivity of the 
individual det{Ma,b), we expect that factor to reduce the evolution of P as compared with the first 
model. Indeed, in the high temperature limit the suppression is through powers of the ratio of the 
evaporative energy scales of the two species, but for low temperatures, det(A^) ~ 2det(M;,). Thus, 
as may have been expected due to the larger thermal inertia of the entire system compared with 
that of an individual species, the actually overall evaporative cooling efficiency is suppressed at low 
temperatures relative to the first model we considered. 



This means that if we cool in the active mode, will be large initially and thus the temperature 
of the entire system will not drop dramatically as one evaporates. Efficient cooling could occur if 
we are able to get to /U ~ /U*, but it is difficult to reach that regime because, although we can 
evaporate all of species "b" (thereby having a continuous curve for the dashed curve in Figure 4 
representing the upper bound for the RHS of Eq.(8)), in the low temperature limit the cooling of the 
whole system proceeds slower, and so the light dotted curve in Figure 4 is roughly twice as steep. 
Thus, the obstruction to reaching the quantum ground state in the two component system in this 
regime (unequal Eevap) can be understood from the considerations and qualitative behavior of the 
first regime (equal Eg^ap)- 

IV: Conclusion: There exists a limit cycle in evaporative cooling a two component fermion system 
that has the consequence of severely limiting the approach to the quantum degenerate ground state. 
One heuristic way to understand this result is that since there is always a gap between the putative 
fermi surface and the evaporative threshold, evaporation can actually "heat up" (that is, disorder) 
the distribution function at the fermi surface. This one way to understand the rather counter- 
intuitive finding that that starting at /i » ji*, the /iP product actually initially decreases during 
evaporation. 

There are many obvious reasons why the two component fermi system is so different than the 
bose system. From the point of view of the discussion here, the trapped bose gas starts with /j, « 
and, analogous to what we have been discussing, evaporatively evolves towards a limit cycle at some 
11* > 0. But, of course due to the singularity of the system at — > 0~, it never gets to //* but 
instead condenses. 

There have been many proposals for surmounting the difficulty of achieving a degenerate fermi 
ground state, and it is not the purpose of this note to review these many inventive ideas. They 
include condensing fermi-bose mixtures (a difficult technical feat) I4,i5,i6,i7,i9 g^^^ using various 
perturbing fields on pure fermi systems. What we suggest is that many cooling techniques that 
evaporate fermions exclusively may be constrained in the same way as described above. 

We would like to end with a brief speculative proposal for reaching lower f3fi products in a 
trapped atomic fermi gas. Consider photoassociating fermi dimers into states just below the trap 
single particle ground state. With fermions there is no stimulated atomic channel back to the trap 
as would the case of photoassociating dimers from a bose condensate; instead pauli blocking and 
the fermi energy both push the system towards dimerization. The remaining fermions then scatter 
off the dimers. In a sense, photoassociating has enhanced the three-body collision rate, which, even 
for identical fermions, is not suppressed by statistics at low energies. Every time a dimer breaks 
in collision, as long as the trapping potential is high enough, the fermions go back into the trap; 
the net effect of creating and breaking dimers in this proposed scheme is to use the difference in 
the dimer pump beams to cool the fermion system "from below" (near the single particle ground 
state of the trap) instead of evaporatively "from above" (that is, above the fermi surface). In that 
sense this scheme has the flavor of bose-fermi mixture schemes, but might be simpler technically. 
Also, this cooling proposal does not a priori require a two fermion mixture, though we imagine 
photoassociating into a dimer composed of dissimilar fermions is likely to be easier than in identical 
ones. It remains to be seen whether such a technique can be practically implemented in a polarized 
atomic fermi system. 

Finally, it would be of great interest to compare the predictions of this simple evaporative model 
directly with experiment. One step in that direction is to generalize the model to include a time 
dependent Eg^ap ^tnd trap lifetime effects. One use of such a direct comparison would be to further 
test how large quantum statistical effects are in current experiments which are far from degeneracy. 
Investigations of such "dynamical" effects are underway but clearly beyond the "kinematic" scope 



and spirit of this note. 
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